Constraining dark energy with Sunyaev-Zel'dovich cluster surveys 
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We discuss the prospects of constraining the properties of a dark energy component, with particular 
reference to a time varying equation of state, using future cluster surveys selected by their Sunyaev- 
Zel'dovich effect. We compute the number of clusters expected for a given set of cosmological 
parameters and propogate the errors expected from a variety of surveys. In the short term they 
will constrain dark energy in conjunction with future observations of type la supernovae, but may 
in time do so in their own right. 



PACS Numbers : 98.80. Es, 98.80. Cq, 98.65.Cw 

Recent observations of type la supernovae (SNe) have 
motivated the search for a ubiquitous energy density 
component, known as dark energy Q. The defining prop- 
erties of this energy are that it has negative pressure and 
does not cluster into galaxies in the same way as dark 
matter, remaining homogeneous on all but the largest 
scales. The standard form is the cosmological constant 
(A), although other possibilities exist including a slowly 
rolling scalar field known as Quintessence. 

The quantification of the properties of this dark en- 
ergy is now a major part of many observational programs. 
One proposal is a satellite, known as SNAP (SuperNova 
Acceleration Probe) |^ which should find around 1800 
SNe out to z « 1.7. This will certainly constrain the 
properties of dark energy |^,^, but without prior infor- 
mation on the matter density, fim, this will have very 
little to say about the time evolution of the equation 
of state parameter = p^/ p^, crucial for distinguish- 
ing between the various dark energy models In this 
Letter, we discuss another approach using future cluster 
surveys selected using the Sunyaev-Zel'dovich (SZ) effect. 
We will show that, dependent on the angular coverage 
(Ai7), frequency {v) and flux limit (6*111,1), such a survey 
may provide complementary information to SNe observa- 
tions, or accurately constrain the properties of the dark 
energy in its own right. 

Observations of clusters via the SZ effect (see ref. 
for a recent review) exploit the fact that the cosmic mi- 
crowave background (CMB) radiation is rescattered by 
hot intracluster gas. Since Compton scattering conserves 
the overall number of photons, the radiation gains energy 
by redistributing them from lower to higher frequencies. 
If one observes them in the Rayleigh-Jcans region of the 
spectrum, the flux of observed photons decreases com- 
pared to the unscattered CMB radiation. The total flux 
depends on the gas mass and mean temperature of the 
cluster, but is independent of their distributions. More- 
over, the number density of such clusters evolves with 
redshift under the action of gravity making it an ideal 
probe of cosmology S. 



The first step is to compute the distribution of clusters 
which will be observed by a particular survey for a given 
set of cosmological parameters. We choose to consider 
the redshift distribution of clusters with mass larger than 
Afiim oc {5m„di/[i/2(l + z)]}3/5i/(z)-2/5 ^hich is given 

by 
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with j^dAI the comoving density of clusters with mass 
between M and M + dM, the volume element, LI{z) 
the Hubble parameter, and dx the angular diameter dis- 
tance. The distribution 4^ could also constrain cosmo- 

dS 

logical parameters and may be powerful if there is sparse 
redshift information available. However, we do not ex- 
pect it to be very sensitive to the equation of state of 
dark energy since the crucial redshift dependence is inte- 
grated out. 

The limiting mass Mjim of the survey can be directly 
related to the total limiting flux Swa^ of the SZ survey by 
the virial theorem and the SZ flux |p|-|ll| . We assume that 
the geometry of the universe is flat and that the late time 
dynamics is dominated by a matter component with den- 
sity rim and a dark energy component with = 1 — rJm. 
Since there are a wide range of dark energy models dis- 
cussed in the literature (see refs. p|,p|jl^] and references 
therein) which all have potentially different late time be- 
haviour we choose to parametrize the equation of state 
by its late time evolution Wrj, = wq + wiz. The comoving 
number density is taken from a series of N-body sim- 
ulations ll^, which yields results similar to using the 
Press-Schechter formalism but predicts more mas- 

sive and less 'typical' clusters, as observed in the simu- 
lations ||l^]. The linear growth factor is computed for a 
given cosmology by solving the ODE for the linear pertur- 
bations numerically and non-linear evolution is taken 
into account via the spherical collapse model. 

In Fig. |l| we illustrate the dependence of the redshift 
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FIG. 1. In the left panel we show the cosmology depen- 
dence of the redshift evolution of the number of clusters and 
in the right panel the mass threshold. We vary f2m, erg, wo 
and wi as explained in the text. The results correspond to 
the experimental setup (II) with 10^ deg^ sky coverage. 



distribution of SZ clusters and the limiting mass on cos- 
mology. The solid line is a model with = 0.3, the 
Hubble constant Hq = 65kmsec~^Mpc~^, as — 0.925, 
Wo = — 1, Wi = and spectral index of density fluctua- 
tions n = 1. The dotted line has fim = 0.5, the dashed 
line CTg = 0.975, the long dashed line wq — —0.8 and the 
dot-dashed line wo — —0.8 and wi — 0.3. The depen- 
dence on n is very weak ||l0| and we therefore fix n = 1. 
We will consider the possible dependence of the number 
density on the parameters 8 = (Hq, as, J^m, wq, wi) in 
the subsequent analysis. From Fig. |l| we see that dN/dz 
is strongly dependent on fim and as, while the depen- 
dence on Wq is still recognizable and that on wi is rela- 
tively weak. 

We make the optimistic assumption that all the clus- 
ters found in the complete surveys can be located suffi- 
ciently well so as to determine their redshift out to some 
critical value Zmax- Furthermore, we will assume that 
this will be known within a precision of Az — 0.01 which 
will allow us to use data bins of size Az. This level of 
accuracy will only require the redshifts to be determined 
photometrically and will be possible using SDSS (Sloan 
Digital Sky Survey) and VISTA (Visible and Infrared 
Survey Telescope for Astronomy) . We can then compare 
to theoretical models using the Cash C statistic [ p^|Jl7| 
for the log-likelihood assuming that the errors are Pois- 
son distributed. 

A number of surveys are expected which are designed 
to detect all clusters above some limiting mass Mii^{z). 
For the purposes of our discussion we will group them 
into four categories whose observing strategies, approx- 
imate Miim and projected number of observed number 
of clusters in a dark energy based cosmology are tabu- 
lated in Table I. The first category (I) of deep and nar- 
row surveys contains the interferometric arrays AMI fl^ ] 
(Arc- minute Micro-Kelvin Imager), SZA ^ (SZ Ar- 
ray) and AMiBA |2^ (Array for Microwave Background 
Anisotropy). For AMI detailed simulations of the sur- 





(I) 


(11) 


(III) 


(IV) 


iS'lini 


0.1 


5 


f»36 




U 


15 


30 


RilOO 




An 


10 


10'' 


20600 


4000 


Miim 


1.5 


^ 7.0 


f» 6.0 


2.5 




fa 90 


^ 1970 


^ 5200 


^ 13600 


SHo 


zbco 


±15 


±15 


-10/ + 5 


5as 


±0.075 


±0.03 


±0.02 


±0.007 




-0.07/+0.10 


±0.05 


±0.03 


±0.02 


Swo 


-OO/+0.30 


-0.15/+0.29 


-0.09/+0.12 


-0.04/+0.12 


5wi 


±CO 


-0.60/-I-0.14 


-0.46/+0.10 


-0.55/+0.05 


SHo 


-15/+5 


±5 


±7 


±5 




-0.04/+0.08 


±0.03 


±0.02 


±0.01 




-0.07/+0.28 


-0.03/+0.14 


-0.09/+0.12 


±0.03 


Swi 


-OO/+0.15 


-0.47/+0.09 


-0.42/+0.06 


±0.03 




±0.04 


±0.02 


±0.01 


±0.01 




-0.07/+0.23 


-0.02/+0.10 


-0.07/+0.04 


±0.03 


5wi 


-OO/+0.15 


-0.40/+0.04 


-0.28/+0.07 


±0.03 



TABLE I. The properties of the different classes of ex- 
periments, the number of clusters one would expect to 
observe in a fiducial cosmology and the 1-cr uncertain- 
ties on the parameters one would deduce for the same 
cosmology if one (a) had no prior information, (b) with 
fixed Sag and (c) imposed both SHq = 5 and fixed 
5as together. The units of (-ffo, i^, An, Mum) are 
(kmsec-^Mpc-\mJy,GHz,deg^ 10"/i"^Mq). We used oo 
to denote cases where we were unable to make a sensible con- 
straint on a particular parameter. 

vey yield have been performed and radio source con- 
tamination has been considered. The second group (II) 
are shallow and wide surveys of which OCRA ||2^ (One- 
Centimetre Receiver Array) is an example. Here, we use 
the flux sensitivity for a single receiver from the proposed 
100 beam array, without taking into account the effects 
of atomospheric water vapour at the site. The third class 
(III), which are shallow but nearly all-sky, correspond to 
what might be possible based on component separation 
using the multi-frequency channels of the PLANCK sur- 
veyor. As an example of the sensitivity we list the 100 
GHz channel. The final category (IV), are deep and wide 
surveys, such as a 1000 element bolometric array which 
may be mounted on a telescope at the south pole. In 
the last case, due to lack of precise figures we will use a 
constant limiting mass [pT| . 

The 1-a errors one would expect on 5Q are 
listed in Table I for a fiducial cosmology 8 — 
(65,0.925,0.3,-0.8,0.3) assuming no prior information 
and 2;,„ax = 1-5. This particular cosmology was cho- 
sen since, firstly, it is in the middle of the parameter 
range preferred by the current data and, secondly, it cor- 
responds to a particular dark energy model which one 
might want to constrain |l^. We have tested the sta- 
bility of our results to small changes in the parameters 
compatible with the current observational data. 



2 






50 





40 


G 







50 





20 




0. 




0.50 
0.40 



G 




0.90 1.00 1.10 



■1.0 -0.5 0.0 



8 



W 







FIG. 2. The marginalized joint likelihood contours in the 
CTg-Slm (left) and wo-f^m (right) planes at the 1-a level. The 
largest contour corresponds to a type (I) survey, the dark grey 
contour to type (II), the light grey contour to type (III) and 
the dashed line contour to type (IV). 

The dependence of ^ on Hq is very weak and the 
double-valued nature of growth factor around w w —0.5 
leads to a degeneracy with the amplitude erg- There- 
fore, it seems sensible to consider the possibility of prior 
assumptions on these two parameters, particularly since 
both should be measured independently of the proper- 
tics of the dark energy by other means. Hq is measured 
using the Hubble Space Telescope at present to within 
AHq — 8 . We will assume that in the next few years 
a precise measurement will be possible to AHq = 5. In 
the case of erg we will assume that it can be measured 
almost exactly by, for example, a low-z X-ray survey. 
Although this will not be precisely true it is useful for 
comparison with ref. The results of imposing the 

prior on erg by itself and combined with that on Hq are 
also listed in Table I. 

There is a clear improvement in one's ability to con- 
strain the parameters in going from a type (I) to type 
(IV). From the point of view of the dark energy the 
salient parameters are flm, wq and wi whose errorbars 
are often asymmetric due to the complicated shape of 
the likelihood surface. Including the prior on ag, appears 
to be useful in removing degeneracies, whereas the dis- 
tribution is very flat in the Hq direction and, therefore, 
inclusion of a prior on it has little significant effect. 

If one uses no prior information with (I), it is only 
possible to measure ug and fl^ accurately and place an 
upper bound on wq. There is no viable constraint on 
wi due to the relatively small number of clusters that 
one would detect in such a setup. If one includes both 
the priors \5Q,ai\ ~ 0.04 and a weak constraint on wq is 
possible, but there is still little information on wi. 

The results of (II) and (III) are qualitatively similar 
with (III) improving on (II). With no prior information 
one can constrain ag and fim considerably (|5(Tg|=0.03, 
|5rJm|=0.05 for (II) and \das\=Om, |^rJm|=0.03 for 
(III)), and good information on wq would be possible. 
However, yet again very little information would be pos- 



sible onwi, a situation which is only mildly alleviated by 
the inclusion of the prior information. It is worth noting 
that for our chosen fiducial model it is easier to set an 
upper bound on wi than a lower bound. This is a general 
trend we observed for the models we studied, though for 
some it was less pronounced. 

Only in the case of (IV) with a fixed value of cg can 
very strong statements be made about wi using this kind 
of observation. Such a setup also gives very accurate 
information on all the other parameters including wq, 
irrespective any prior. This provides clear motivation for 
considering the feasibility of this setup. 




FIG. 3. The 1-cr joint likelihoods in the wq-wi plane. In the 
left panel for setup (II) where dark shaded region is obtained 
with a maximum redshift of Zmax = 1.5 and the light shaded 
region corresponds to a maximum redshift of Zmax = 0.5. In 
the right panel for setups (II), (III) and (IV) using the same 
conventions as in Fig. ^. The transparent 3-dot dashed line 
corresponds to the joint likelihood for the SNe survey SNAP. 

We have already noted that the errorbars are in gen- 
eral very asymmetric. In order to investigate this we have 
plotted in Fig. ||the joint likelihood surfaces in the crg-fim 
and wo-rim planes for each of the setups (I)- (IV), which 
show visually the relative improvement that one might 
expect. The degeneracies are similar to those observed 
in previous work and we see that only (II) , (HI) and 
(IV) constrain wq in any significant way. Nonetheless, it 
is clear that in each case the value of fim is constrained 
extremely well. We have used ^max = 1.5, however, the 
using Zjiiax = 0.5 remarkably has little effect on the size 
of the errorbars, since it is the statistical weight of the 
large number of clusters found at low rcdshifts which 
fixes these parameters. We also performed an analysis 
with Az ~ 0.025 and found that this increases the un- 
certainties on the estimated parameters in a similar way 
to changing z^ax = 1.5 to z,nax = 0.5. 

The degeneracy between wq and wi is particularly im- 
portant from the point of view of dark energy and this 
is illustrated in Fig. ^, left panel, for (II) under the opti- 
mistic assumption that Zmax = 1.5 and when z^ax = 0.5. 
The degeneracy has a complicated, double-valued shape 
and the constrained region is much smaller when Zmax 
is larger. This is as expected since constraining these 
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parameters requires more information at high redshifts. 

Our results show that only for setup (IV) and an ef- 
fectively fixed value of erg can one independently fix the 
crucial parameter wi using this kind of measurement. 
However, all is not lost; it was pointed out in ref. Q that 
given independent prior information on Q^, SNe mea- 
surements can accurately constrain the dark energy. As 
we have already pointed out even setup (I) will provide 
important information in this respect and the others will 
improve on this. 

Even more information can be gleaned by making the 
comparison of the two different probes of dark energy in 
the wq-wi plane. Fig. |[ right panel, illustrates this for 
setups (II), (III) and (IV) compared to a similar calcu- 
lation for SNAP Even for (II) performances of the 
two methods are comparable in terms of the area of the 
1-a contour and for (IV) the result is very much better. 
Notice also that the degeneracy in this plane is also to- 
tally different and combining them would give a localized 
region pinning down wq very accurately and wi to within 
^ ±0.2. While this may not be enough to rule out wi = 
at the 2-a level, a look at the various models for dark en- 
ergy considered in ref. shows that such observations 
would put tight constraints on the particular dark energy 
models. 

Our basic philosophy has been to investigate the abso- 
lute best case constraints that a given survey can achieve 
in terms of the properties of the dark energy. In this 
spirit, we have shown that as with SNe observations, 
cluster surveys selected using the SZ effect will provide 
important information as to the nature of the dark en- 
ergy and that there is a potential synergy between the 
two. However, our conclusions were drawn from a highly 
idealized model of cluster physics. 

One of the key sources of systematic uncertainties will 
come from the — 5*11,11 relation due to, for exam- 

ple, heat input or the clusters being not completely viri- 
alized |^^. These effects might manifest themselves in 
terms of either a systematic shift in the overall normaliza- 
tion, or in statistical scatter. We have estimated the pos- 
sible effects of these uncertainties ||ll|] and found that if 
the scatter is less than about 20% and the overall normal- 
ization is accurate to with 5% then one can distinguish 
our fiducial model from the standard ACDM model. To 
get an idea of why the constraint on required accuracy 
of the overall normalization is particularly important we 
just comment that a 20% change would lead to a factor 
of 2 change in the total number of clusters. 

We are optimistic that many of the practical difficulties 
which we have ignored will be addressed with the first 
generation of SZ survey instruments and can be taken 
into account in the future with the qualitative picture of 
our results will remain: that SZ cluster surveys provide 
a robust complementary probe for dark energy. 
ACKNOWLEDGMENTS: We thank R. Crittenden, 
M. Jones and P. Wilkinson for useful discussions. The 



parallel computations were done at the UK National 
Cosmology Supercomputer Center funded by PPARC, 
HEFCE and Silicon Graphics / Cray Research. JW and 
RB are supported by a PPARC. 



[1] S. Perlmutter et al, Ap. J. 483, 565 (1997). A. Riess et 
al, Astron. J. 116, 1009 (1998). S. Perlmutter et al, Ap. 
J. 517, 565 (1999). 

[2] C. Wetterich, Nucl. Phys. B302, 668 (1988). B. Ratra 
and P. J. E. Peebles, Phys. Rev. D37, 3406 (1988). 

[3] snap.lbl.gov 

[4] D. Huterer and M. S. Turner, Phys. Rev. D 60, 081301 
(1999). I. Maor, R. Brustein and P. J. Steinhardt (2001), 
Phys. Rev. Lett. 86, 6. Erratum: Phys. Rev. Lett. 87, 
049901. P. Astier, Phys. Lett. B 500, 8 (2001). 

[5] J. Weller and A. Albrecht, Phys. Rev. Lett. 86, 1939 
(2001). J. Weller and A. Albrecht, astro-ph/0106079. 

[6] R. A. Sunyaev and Ya. Zel'dovich, Comm. Astro- 
phys. Space Phys. 4, 173 (1972). R. A. Sunyaev and 
Ya. Zel'dovich, MNRAS 190, 143 (1980). 

[7] M. Birkinshaw, Phys. Rep. 310, 97 (1999). 

[8] P. B. Lilje, Ap. J. Lett. 386, L33 (1992). J. Oukbir and 
A. Blanchard, Astron. and Astro. 262, L210 (1992). V. 
Eke, S. Cole and C. Frenk, MNRAS 282, 263 (1996). 
N. A. BahcaU and X. Fan, Ap. J. 504, 1 (1998). 

[9] P. T. P. Viana and A. R. Liddle, MNRAS 281, 323 
(1996). L. Wang and P. J. Steinhardt, Ap. J. 508, 483 
(1998). 

[10] Z. Haiman et al, Ap. J. 553, 545 (2001). 
[11] J. Weller, R. A. Battye and R. Kneissl, In preparation. 
[12] Ph. Brax and J. Martin, Phys. Lett. B468, 40 (1999). 
[13] A. Jenkins et al, MNRAS 321, 372 (2001). 
[14] W. H. Press and P. Schechter, Ap. J. 187, 425 (1974). 
[15] E. Pierpaoh, D. Scott and M. White, MNRAS 325, 77 
(2001). 

[16] W. Cash, Ap. J. 228, 939 (1979). 

[17] G. Holder, Z. Haiman and J. Mohr, Ap. J. Lett. 560, 111 
(2001). 

[18] R. Kneissl et al, MNRAS 328, 783 (2001). M. E. Jones, 
astro-ph/0109351. 

[19] G. P. Holder et al., Ap. J. 544, 629 (2000). J. J. Mohr et 
al., in 'Extrasolar Planets to Cosmology: The VLT Open- 
ing Symposium' (2000), edited by A. Renzini, (2000). 
J. E. Carlstrom et al., in 'Constructing the Universe with 
Clusters of Calaxies' (2000), edited by F. Durret and 
G. Gerbal, (2000). 

[20] K. Y. Lo et al., in 'New Cosmological Data and the Val- 
ues of the Fundamental Parameters' (2000), edited by 
A. Lasenby and A. Wilkinson, (2000). 

[21] I. W. Browne et al, in 'Radio Telescopes' (SPIE Proc. 
Vol 4015), edited by H. R. Butcher, (2000). 

[22] W. L. Freedman et al., Ap. J. 553, 47 (2001). 

[23] L. Verde, Z. Haiman and D. N. Spergel, astro- 
ph/0106315. 



4 



